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Abstract 

We report the first experimental study of the quantum interference correction to the conductivity of bilayer graphene. 
Low-field, positive magnetoconductivity due to the weak localisation effect is investigated at different carrier densities, 
including those around the electroneutrality region. Unlike conventional 2D systems, weak localisation in bilayer graphene 
is affected by elastic scattering processes such as intervalley scattering. Analysis of the dephasing determined from the 
magnetoconductivity is complemented by a study of the field- and density-dependent fluctuations of the conductance. Good 
agreement in the value of the coherence length is found between these two studies. 
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1. Introduction 

The recently developed technique to fabricate free- 
standing sheets of graphene [1] brought with it the first 
opportunity to study a true two-dimensional crystal. 
Since this development, remarkable phenomena have 
been predicted and observed in both monolayer and 
bilayer forms of this crystal [2]. One of the important 
questions that has arisen is the influence of the chiral 
nature of charge carriers in graphene on the quantum 
interference. In recent theories [3,4], this influence was 
considered on the quantum correction to the conduc- 
tivity due to weak localisation ( WL) . This phenomenon 
is well studied in conventional 2D systems [5], but in 
graphene the quantum interference was predicted to 
be affected not only by inelastic scattering, but also by 
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elastic inter- and intra- valley scattering. It was shown 
that the shape of the low-field magnetoconductivity 
(MC) of bilayer graphene [4] can reflect the interplay 
between these scattering mechanisms. 

Here we report an experimental investigation of the 
MC of bilayer graphene. We show how the different 
scattering processes are revealed in the weak localisa- 
tion and how they depend on the carrier density. 



2. Experiment 

The sample was fabricated by the procedure devel- 
oped in [1] . Highly-oriented pyrolytic graphite was me- 
chanically exfoliated and deposited on top of a ther- 
mally oxidized n + Si substrate (with an oxide thick- 
ness of 300 nm). A pair of Ohmic (Au/Cr) contacts 
was deposited using electron-beam lithography on op- 
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posite edges of a bilayer flake of length 1.5 /im and 
width 1.8 fim. The highly-doped substrate was used as 
a back gate and the carrier density in the flake was 
controlled by the gate voltage V g . (Using the formula 
for a plane-plate capacitor, one can find a simple re- 
lation between the gate voltage and the carrier den- 
sity: dn/dVg = 7.3 x 10 10 cm" 2 V _1 .) Electrical mea- 
surements were performed in the temperature range 
from 0.04 K to 4.2 K. The gate voltage dependence of 
the resistance of the studied sample is presented in 
Fig. 1(a). Far from electroneutrality region at V g = 
0, where the charge carriers change from electrons to 
holes, the sample exhibits an almost constant mobility 
of » 7000cm 2 /V-s. 

To verify that we are dealing with a true bilayer 
graphene system the quantum Hall effect was studied. 
The results are shown in Fig. 1(b). The values of the 
conductance at the plateaux correspond to those of 
a conventional 2D system with four-fold degeneracy; 
however, the filling factor v at which the first plateau 
occurs is unconventional, and is a result of the presence 
of the lowest Landau level at zero energy. This is a clear 
demonstration that the sample is bilayer graphene [2] . 
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Fig. 1. (a) Field-effect in bilayer graphene. The numbered 
boxes indicate the regions where the magnetoconductivity 
was studied, (b) Quantum Hall effect in bilayer graphene. 
(c), (d) Temperature smearing of conductance fluctuations 
as a function of gate voltage (carrier density). 

The low-field MC of the sample is strongly affected 
by fluctuations of the conductance, seen in Fig. 1. Since 
the size of the sample is rather small, these repro- 
ducible, mesoscopic fluctuations are significant. We ob- 
serve fluctuations with amplitude ~ e 2 /h, both as a 
function of the gate voltage (Fig. l(c,d)) and magnetic 
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Fig. 2. Conductivity of the sample at low magnetic field 
near electroneutrality region (region I, top panel) and at 
higher concentration (region II, bottom panel). 

field (Fig. 2). This amplitude decreases with increas- 
ing temperature, Fig. l(c,d), which is comparable to 
the decrease of the amplitude of the MC. Therefore, 
straightforward analysis of the raw MC data is not pos- 
sible and we have to apply an averaging procedure. We 
averaged the effect of the magnetic field on the conduc- 
tivity across each of the studied regions, highlighted by 
boxes in Fig. 1(a) (AV g = 2 V). (A detailed descrip- 
tion of this averaging procedure is given in [6].) Results 
of the averaging are presented in Fig. 3 where the low- 
field MC in the electroneutrality region and electron 
region is shown. The experimental points are fitted by 
the theoretical expression [4]: 
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Here F(z) = In z + ^/>(0.5 + z~ ) , ip(x) is the digamma 
function, Tg 1 = AeDB/h, t^ 1 is the phase-breaking 
rate and r^T 1 = t^ 1 + r^ 1 , where t w is the intra- valley 
"warping" time and t z is the time of chirality breaking 
[7] . It is seen from Eq. 1 that the magnetoconductivity 
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is affected not only by inelastic scattering (r^ ) but 
also by elastic scattering rates r" 1 and t^ 1 . 
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Fig. 3. Low-field magnetoconductivity of bilayer graphene. 
Solid lines are fits to the data by Eq. 1. The dashed line is a 
fit by the same equation if only the first term is considered. 

The solid lines in Fig. 3 present the best fit of the ex- 
perimental points by the first and second terms of Eq. f . 
The third term is found to be significantly smaller than 
the other two and produces only a minor correction to 
the fit. The second term plays a crucial role: the satura- 
tion of the WL correction at higher fields is due solely 
to this term. If only the first term is considered, the MC 
increases continuously with increasing magnetic field. 
This is shown by the dashed line in Fig. 3, which is a 
best fit to the low-field data by the first term alone. 



3. Discussion 



than in the electron region. A possible reason for this is 
that in the electroneutrality region the carrier density, 
which is on average zero, becomes inhomogeneous due 
to the formation of electron-hole puddles. This results 
in an effective reduction of the size of the sample. 

We analyse the fluctuations of the conductance in 
terms of the conventional theory of universal conduc- 
tance fluctuations [8]. From the amplitude of these fluc- 
tuations, averaged over the magnetic field range AB = 
0.5 T, we determine the coherence length as a function 
of the temperature, Fig 4(b). These values of Z/0 agree 
well with the results of the WL analysis. The elastic 
scattering in bilayer graphene therefore appears to af- 
fect the fluctuations of the conductance less than the 
magnetoconductivity due to WL. 

Figure 4(c) shows the relation between the char- 
acteristic lengths in bilayer graphene: the saturation 
value of 1/0, the intervalley scattering length Li, and 
the mean free path l p . As mentioned above, the satu- 
ration value of 1/0 is seen to be somewhat smaller in 
the electroneutrality region than at higher densities, 
as has also been seen in monolayer graphene systems 
[9]. In all studied regions, Li is found to be temper- 
ature independent and comparable to the saturation 
value of Z/0. This implies that, similar to [9], there is 
a significant contribution of the sample edges to in- 
tervalley scattering. There is a trend of a decrease of 
Li with increasing density; but its origin is not clear, 
as, contrary to single-layer, the density of states in bi- 
layer graphene is energy independent. Both L^ and Li 
are found to be significantly larger than the mean free 
path l p , which justifies the application of the diffusive 
theory [4] . The small contribution of the third term in 
Eq. 1 is in agreement with theoretical expectations [4]: 
T w ~ t p for the studied range of densities. 

The temperature dependence of the dephasing rate is 
presented in Fig. 4(d). At temperatures > 1 K it can be 
described by the conventional theory [5] for dephasing 
due to the electron-electron interaction in the diffusive 
regime, Tr p < 1: 



Figure 4 contains the main results of our analysis. 
Figure 4(a) presents the temperature dependence of 
the coherence length L^ for the electroneutrality and 
electron regions. At low temperatures, 1/0 saturates at 
a value comparable to the size of the sample. This sat- 
urated value is smaller in the electroneutrality region 
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where g — a/(e 2 /h) is a dimensionless conductivity, 
and P is an empirical coefficient close to unity. 
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Fig. 4. (a) Coherence length found from the WL analysis 
for regions I and II. (b) Values of found from the anal- 
ysis of the amplitude of the conductance fluctuations as a 
function of the magnetic field, (c) Comparison of the char- 
acteristic lengths in bilayer graphene at different concen- 
trations: saturation value of (open circles), intervalley 
scattering length Li (half-filled circles) and mean free path 
l p (filled circles). (The dashed line is is a guide to the eye 
to show the decrease of the saturation value of in the 
electroneutrality region.) (d) Temperature dependence of 
the dephasing rates in the same regions. Lines are drawn 
according to Eq. 2. 
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4. Conclusion 



The quantum correction to the conductivity due to 
weak localisation in bilayer graphene is affected both 
by inelastic and elastic scattering mechanisms. This 
is caused by the chirality of charge carriers in this 
system. Values of the coherence length and interval- 
ley scattering length are determined from analysis of 
the WL. The found values of the coherence length are 
in agreement with those determined from analysis of 
the observed conductance fluctuations as a function of 
magnetic field. The latter is done in terms of the con- 
ventional, unmodified theory of universal conductance 
fluctuations. The temperature dependence of the de- 
phasing rate at high temperatures can be described by 
the theory of electron-electron interaction in the dif- 
fusive regime. 
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